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Abstract. We provide tables of theoretical isochrones in several photometric systems. To this aim, the following 
steps are followed: (1) First, we re-write the formalism for converting synthetic stellar spectra into tables of 
bolometric corrections. The resulting formulas can be applied to any photometric system, provided that the zero- 
points are specified by means of either ABmag, STmag, VEGAmag, or a standard star system that includes 
well-known spectrophotometric standards. Interstellar absorption can be considered in a self-consistent way. (2) 
We assemble an extended and updated library of stellar intrinsic spectra. It is mostly based on "non-overshooting" 
ATLAS9 models, suitably extended to both low and high effective temperatures. This offers an excellent coverage 
of the parameter space of T e fj, log g, and [M/H] . We briefly discuss the main uncertainties and points still deserving 
more improvement. (3) From the spectral library, we derive tables of bolometric corrections for Johnson-Cousins- 
Glass, HST/WFPC2, HST/NICMOS, Washington, and ESO Imaging Survey systems (this latter consisting on 
the WFI, EMMI, and SOFI filter sets). (4) These tables are used to convert several sets of Padova isochrones into 
the corresponding absolute magnitudes and colours, thus providing a useful database for several astrophysical 
applications. All data files are made available in electronic form. 
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1. Introduction 

One of the primary aims of stellar evolution theory is 
that of explaining the photometric data - eg. colour- 
magnitude diagrams (CMDs), luminosity functions, colour 
histograms - of resolved stellar populations. To allow 
whatever comparison between theory and data, the basic 
output of stellar models - the surface luminosity L and ef- 
fective temperature T e s - must be first converted into the 
observable quantities, i.e. magnitudes and colours. This 
conversion is performed by means of bolometric correc- 
tions (BC) and T e ff-colour relations, and later by consid- 
ering the proper distance, absorption and reddening of the 
observed population, and the photometric errors. 
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Determining BCs and T e ff-colour relations is indeed 
one of the most basic tasks in stellar astrophysics. 
Empirical determinations (see e.g. the compilations by 
Schmidt-Kaler 1982, Flower 1996, and Alonso et al. 
1999ab) involve a great observational effort, and are obvi- 
ously the most reliable if compared to purely theoretical 
determinations. However, since the empirical calibrations 
are mostly based on nearby stars, only a limited region in 
the space of stellar parameters (T e g, log g, and metallic- 
ity [M/H]) can be covered by observations. This contrasts 
with our present-day capabilities of getting resolved pho- 
tometry of populations that are certainly very different 
from the local one, like those in dwarf galaxies and in the 
Bulge. For instance, present empirical relations do not in- 
clude young- metal poor populations, or super-metal rich 
stars, which are likely present in the resolved galaxies of 
the Local Group. 
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Another important limitation of empirical relations is 
that they are usually available just for a small set of fil- 
ters or photometric systems, like the popular Johnson- 
Cousins-Glass one. Again, this is in contrast with the rapid 
diffusion in the use of specific filter sets, for which no em- 
pirical relation is yet available, though large databases are 
already being collected. Just to mention a few relevant ex- 
amples, new filter sets have been adopted in the Hubble 
Space Telescope (HST) Wide Field Planetary Camera 2 
(WFPC2), in the European Southern Observatory (ESO) 
Wide Field Imager (WFI), and in the Hipparcos mission. 
Moreover, brand-new photometric systems have been de- 
signed for the Sloan Digital Sky Survey (SDSS), and will 
be adopted in the future GAIA mission. An impressive 
amount of data will be provided by these instruments in 
the coming years, and much of it will soon become of pub- 
lic access. Then, it would be highly desirable to have the 
capability of converting stellar models to these many new 
photometric systems, bypassing for a moment the time- 
consuming procedure of empirical calibration. 

With this target in mind, we have undertaken a 
project aiming at providing theoretical BCs and colour 
transformations for any broad-band photometric system, 
and many intermediate-band systems as well. Actually, 
the project starts with the work by Bertelli et al. 
(1994), who presented a large database of theoretical 
isochrones and converted them to the Johnson-Cousins- 
Glass UBVRIJHK system. The transformations were 
primarily based on Kurucz (1993) synthetic atmosphere 
models, suitably extended in the intervals of lower and 
higher effective temperatures. In Chiosi et al. (1997), the 
same theoretical isochrones are converted to the WFPC2 
photometric system, and the attention is paid to the fea- 
tures that isochrones and single stellar populations present 
in ultraviolet (UV) colours. In particular, they address the 
question whether the variation of UV colours of elliptical 
galaxies as a function of red-shift presents signatures from 
which one can infer the age and type of the source emit- 
ting the UV flux. Later on, Salasnich et al. (2000) present 
new isochrones for a-enhanced chemical mixtures, for both 
UBVRIJHK and WFPC2 photometric systems. 

The present work is a natural follow-up of this project, 
in which, besides extending the number of available pho- 
tometric systems, we aim at updating and improving the 
database of stellar spectra which is at the basis of the 
complete procedure. 

The plan of this paper is as follows: Sect. || details 
the adopted formalism. In Sect. ^, we describe the stellar 
spectral library in use. Sect. | gives the basic informa- 
tion about the several photometric systems under con- 
sideration. The resulting tables of bolometric corrections 
are then applied to a large database of stellar isochrones, 
which are already described in published papers and 
briefly recalled in Sect. H. Sect. |(] illustrates some main 
properties of the derived isochrones, and describes the re- 
trieval of data in electronic form. 



2. Synthetic photometry 

By synthetic photometry we mean the derivation of pho- 
tometric quantities based on stellar intrinsic (and mostly 
theoretical) spectra, rather than on actual observations. 
The first works in this field (e.g. Buser & Kurucz 1978; 
Edvardsson & Bell 1989) were based on sets of synthetic 
spectra covering very modest - by present standards - in- 
tervals of effective temperature, gravity, and metallicity. 
The situation dramatically improved with the release of 
a large database of ATLAS9 synthetic spectra by Kurucz 
(1993). The first systematic use of Kurucz spectra on the- 
oretical isochrones has been from Bertelli et al. (1994) and 
Chiosi et al. (1997). 

In order to apply synthetic photometry to sets of the- 
oretical isochrones, the basic step consists in the deriva- 
tion of bolometric corrections and temperature-colour re- 
lations from the available spectra. Several papers deal with 
the problem (e.g. Bertelli et al. 1994; Chiosi et al. 1997; 
Lejeune et al. 1997; Bessell et al. 1998), presenting math- 
ematical formalisms that, although looking somewhat dif- 
ferent, should be equivalent and produce the same results 
when applied to the same sets of spectra, filters, and zero- 
points. In the following, we re-write this formalism in a 
very simple way. Our aim is to have generic formulas that, 
by just minimally changing their input quantities, can be 
applied to a wide variety of photometric systems. 



2.1. Basic concepts 

For a star, the spectral flux as it arrives at the Earth, fx, 
is simply related to the flux at the stellar surface, Fx, by 



fx 



W-° AA *{R/d) 2 F x 



(1) 



where R is the stellar radius, d is its distance, and Ax is 
the extinction in magnitudes at the wavelength A. Once fx 
is known, the apparent magnitude ms x , in a given pass- 
band with transmission curve S\ comprised in the interval 
[Ai, A 2 ], is given by 



m Sx = -2.5 log 



Sx* ^f\Sxd\ s 
J^XfxSxdX, 



(2) 



where /° represents a reference spectrum (not necessar- 
ily a stellar one) that produces a known apparent mag- 
nitude m? s . In other words, /° and mfg completely de- 
fine the zero-points of a synthetic photometric system (see 
Sect. [D| below). 

In Eq. (||), the integrands XfxSx are proportional to 
the photon flux (i.e. number of photons by unit time, 
surface, and wavelength interval) at the telescope detec- 
tor. This kind of integration applies well to the case of 
modern photometric systems that have been defined and 
calibrated using photon-counting devices such as CCDs 
and IR arrays. However, more traditional systems like 
the Johnson-Cousins-Glass U BV RI J H K LM N one, have 
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been defined using energy-amplifier devices. In this latter 
case, energy integration, i.e. 



-2.5 log 



X 2 f°xSxdX / 



(3) 



would be more appropriate to recover the original sys- 
tem. Notice, however, that the difference between energy 
and photon integration is usually very small, unless the 
pass-bands are extremely wide. Unless otherwise stated, 
in papers of this series we will adopt the integration of 
photon counts. 

2.2. Deriving bolometric corrections 

The starting point of our work are extended libraries of 
stellar intrinsic spectra F\, as derived from atmosphere 
calculations for a grid of effective temperatures T e ff, sur- 
face gravities g, and metallicities [M/H]. The particular 
library we adopt will be described in Sect. |3|. 

From this library, we aim to derive the absolute mag- 
nitudes Ms A for each star of known (T e s,g, [M/H]) - and 
hence known Fx ■ This can be obtained by means of Eq. (||) 
(or 0), once a distance of d = 10 pc is assumed, i.e. 



M 



Sx 



-2.5 log 



R 



10 pc 



: f£ \FxlO-° AA *Sxd\ 



dX 



4(4) 



and once the stellar radius R is known. Since the quantities 
(T c ff, <7, [M/H]) are not enough to specify R (in this case 
we need to have also the stellar mass M), we should first 
eliminate R from our equations. This is possible if we deal 
with the bolometric corrections, 



BC, 



Sx 



Mi 



bol 



Mc 



(5) 



From the definition of bolometric magnitude, we have (see 
also Bessell et al. 1998 for a similar approach): 



Mboi = M bo i, - 2.5 log(L/L Q ) 

= M bol , - 2.5 log(4 7 r J R 2 J F bo i/ J L ) 



(6) 



where F bo i = J °° Fxd\ = (tT^ s is the total emerging flux 
at the stellar surface. 

Substituting Eqs. (§) and (@) into Eq. (§), we get 



BC S 



M bol , - 2.5 log [4tt(10 pc) 2 F bol /L Q ] 
'f^ 2 \F X 10-°- 4A ^ A dA\ 



(7) 



+2.5 log 



that, as expected, depends only the spectral shape 
(F\ 10 _0 - 4 ' 4A /Fboi), and on basic astrophysical constants. 

To keep consistency with our previous works (e.g. 
Salasnich et al. 2000), we adopt M bo i,© = 4.77, and 
L = 3.844 x 10 33 ergs -1 (Bahcall et al. 1995). 

By means of Eq. (0) , we tabulate BCs x for all spectra 
in our input library, and for several different photometric 
systems. The BC$ X can be then derived for any intermedi- 
ate (Tctf^g, [M/H]) value, by interpolation in the existing 



grid. We adopt simple linear interpolations, with logT e ff, 
logg, and [M/H] as the independent variables. 

Next, to attribute absolute magnitudes to stars of 
given (T c ff,L) along an isochrone, we simply compute 
Mboi with Eq. (||), and hence 

Msx = M bo i - BC Sx • (8) 

In this formalism, an extinction curve Ax can be ap- 
plied to all spectra of the stellar library, so as to allow the 
derivation of bolometric corrections (and synthetic abso- 
lute magnitudes) that already include extinction in a self- 
consistent way (cf. Eqs. [?] and||). With this approach, the 
extinction on each pass-band depends not only on the to- 
tal amount of extinction, but also on the spectral energy 
distribution of each star - i.e. on its spectral type, lumi- 
nosity class, and metallicity (see also Grebel & Roberts 
1995). Anyway, for the sake of simplicity, in the present 
work we will deal with the case Ax = only. Tables for dif- 
ferent extinction curves will be discussed in a forthcoming 
paper. 

Finally, it is worth mentioning that the newly-defined 
SDSS photometric system makes use of an unusual defini- 
tion for magnitudes (see Lupton et al. 1999). This specific 
case, for which some of the above equations do not apply, 
will also be discussed in a subsequent paper of this series. 

2.3. Reference spectra and zero-points 

By photometric zero-points, one usually means the con- 
stant quantities that one should add to instrumental mag- 
nitudes in order to transform them to standard magni- 
tudes, for each filter Sx- In the formalism here adopted, 
however, we do not make use of the concept of instrumen- 
tal magnitude, and hence such constants do not need to be 
defined. Throughout this work, instead, by "zero-points" 
we refer to the quantities in Eqs. (||) and (||) that depend 
only on the choice of /° and m s . They are constant for 
each filter, and are responsible for the conversion of the 
synthetic magnitude scale into a standard system. 

As for these quantities, there are four different cases 
of interest: 

2.3.1. VEGAmag systems 

They make use of Vega (aLyr) as the primary calibrat- 
ing star. The most famous among these systems is the 
Johnson-Cousins-Glass UBVRIJHKLMN one, that can 
be accurately recovered by simply assuming that Vega has 
V = 0.03 mag, and all colours equal to 0. Other systems, 
like Washington and the HST/WFPC2 VEGAmag one, 
follow a similar definition (some colours, however, are de- 
fined to have values slightly different from 0). 

Calibrated empirical spectra of Vega are available (e.g. 
Hayes & Latham 1975; Hayes 1985), covering the wave- 
length range from 3300 to 10500 A, an interval that can 
be extended up to 1150 A when complemented with IUE 
spectra (Bohlin et al. 1990). They can be used to de- 
fine VEGAmag systems in the optical and ultraviolet. 
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However, as the wavelength range accessible to present 
instrumentation is much wider, a Vega spectrum cover- 
ing the complete spectral range has become necessary. 
Synthetic spectra as those computed by Kurucz (1993) 
and Castelli & Kurucz (1994), fulfill this aim. In this case, 
the predicted fluxes at Vega's surface, Fa, are scaled by 
the geometric dilution factor 



(R/d) 2 = (0.5 6» d /206264.81) 2 



(9) 



where 9d is the observed Vega's angular diameter (in arc- 
sec) corrected by limb darkening. 

More recently, composite spectra of Vega have been 
constructed by assembling empirical and synthetic spec- 
tra together (e.g. Colina et al. 1996), so that some small 
deficiencies characteristic of synthetic spectra^] are cor- 
rected. This has the precise scope of providing a reference 
spectrum for conversions between apparent magnitudes of 
real (observed) stars, and physical fluxes. However, it is 
not clear whether such composite Vega spectrum should 
be preferable when synthetic photometry is performed on 
theoretical spectra, as in the present work. For instance, 
if the ATLAS9 spectrum for Vega has the core of Bahner 
lines differing by as much as 10 percent from the ob- 
served ones (cf. Colina et al. 1996), it is probable that 
the same deviations will be present in all Kurucz spectra 
of comparable temperatures/gravities. If this is the case, 
the synthetic Vega spectrum would probably give better 
zero-points for these stars than the composite Colina et 
al. (1996) one. 

For this reason, we simply adopt the synthetic 
ATLAS9 model for Vega, with T eS = 9550 K, log 5 = 3.95, 
[M/H] = —0.5, and microturbulent velocity £ — 2kms _1 , 
the spectrum being provided by Kurucz (1993). Castelli 
& Kurucz (1994) have computed a higher-resolution 
spectrum for the same model, using the more refined 
ATLAS 12 code. As discussed by these authors, ATLAS9 
and ATLAS 12 spectra for Vega are almost identical. 

Once the synthetic model i^ ega is chosen, we just need 
to adopt a fixed value for the dilution factor (R/d) 2 , in 
order to have / A ° ga at the Earth's surface (outside the 
atmosphere). Two choices are possible then: We either (i) 
adopt an observed value of 9d as input to Eq. (||) , or (ii) 
adopt the observed Vega flux as measured at the stellar 
surface, at a given wavelength, as an input to Eq. ([!]). 

Since direct measures of 9d are relatively more uncer- 
tain than direct measures of f^ es& , we prefer to adopt the 
second alternative: Taking the flux values at 5556 A from 
Hayes (1985; /J oga = 3.542 x lO" 20 ergs" 1 cm" 2 Hz" 1 ), 
and at 5550 A from Kurucz (1993) Vega model (F^ cga = 
5.507 x 10 7 ergs" 1 cm" 2 A" 1 ), we obtain (R/d) 2 = 6.247 x 
10~ 17 . This value implies an angular diameter of 9d = 
3.26 mas (Eq. ^) for Vega, that compares very well with 
the observed values of 3.24 ± 0.07 mas (Code et al. 1976) 
and 3.28 ± 0.06 mas (Ciardi et al. 2000). 

For a discussion of the differences between synthetic and 
observed Vega spectra, the reader is referred to Castelli & 
Kurucz (1994) and Colina et al. (1996). 



It is worth remarking that, since the zero-points in 
VEGAmag systems are attached to the observed Vega 
fluxes, their synthetic absolute magnitudes may have sys- 
tematic errors of a few hundredths of magnitude (say up 
to 0.03 mag), which is the typical magnitude of errors in 
measuring fluxes at the Earth's surface. Somewhat smaller 
errors, however, are expected in the colours. These uncer- 
tainties will probably not be eliminated unless definitive 
9d and fY° sa measurements become available. 



2.3.2. ABmag systems 

In the original work by Oke (1964), monochromatic AB 
magnitudes are defined by 



toab,^ = -2.5 \ogf u -48.60 . 



(10) 



This means that a reference spectrum of constant flux 
density per unit frequency 



JA 



AB,i/ 



3.631 x 10 -20 erg s" 1 cm" 2 Hz" 



(11) 



will have AB magnitudes ™ AB v ~ a ^ frequencies v. 

This definition can be extended to any filter system, 
provided that we replace the monochromatic flux /„ with 
the photon counts over each pass-band S\ obtained from 
the star, compared to the photon counts that one would 
get by observing / ABl/ : 



™ab,s a 



-2.5 log 



J^(X/hc)f x S x dX • 
ji 2 (A/M/ A B,A<SAdA 



(12) 



where / AB A = / AB V c/X 2 . Then, it is easy to show that 
Eq. ( |l2| ) is is just a particular case of our former Eq. (||), 
for which / AB v is the reference spectrum, and m AB Sx = 
are the reference magnitudes. 

It is worth mentioning that different practical imple- 
mentations of the ABmag system have been defined over 
the years (e.g. Oke & Gunn 1979, and Fukugita et al. 
1996). They differ only in the definition of the reference 
stars (or reference stellar spectra) used as spectrophoto- 
metric standards during the conversion from observed in- 
strumental magnitudes into fluxes f v (that are used in 
Eq. |ici| ) . Since we are dealing with synthetic spectra only, 
such a conversion is not necessary in our case. It follows 
that these different definitions are not a point of concern 
to us. 



2.3.3. STmag systems 

ST monochromatic magnitudes have been introduced by 
the HST team, and are defined by 



tost,a = -2.5 log f x -21.10 



(13) 



This means that a reference spectrum of constant flux 
density per unit wavelength 



/° T A = 3.631 x 10~ 9 erg s" 1 cm" 2 A" 1 



(14) 
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will have ST magnitudes nig T A = at all wavelengths. 
Similarly to the case of AB magnitudes, this can be gen- 
eralized to any pass-band system with: 



J^(X/hc)hS x dX 
f£(\/hc)fg T>x S x d\ 



(15) 



Again, this situation is easily reproduced in our for- 
malism with the adoption of /° = /g T A and Wg T Sx = 0. 



2.3.4. "Standard stars" systems 

This class comprises all photometric systems that use 
standard stars different from Vega to define the zero- 
points. Good examples are the Vilnius (Straizys & 
Zdanavicius 1965) and Thuan-Gunn (Thuan & Gunn 
1976) systems. 

In this case, we are forced to use empirical spectra of 
standard stars in order to define /° and m^ A . A good set is 
provided by the four metal-poor subdwarfs BD +17°4708, 
BD+26°2606, HD 19445 and HD 84937, which are widely- 
used spectrophometric secondary standards (Oke & Gunn 
1983), as well as standards for several photometric sys- 
tems. 

Actually, the present work does not deal with any 
"standard stars" system, and this case is here included 
just for the sake of completeness. Details about a few spe- 
cific systems - including the possible choices for /° - will 
be given in forthcoming papers. 



3. The stellar spectral library 

In the following, we will present the stellar spectral librar 
put together for this work. For the sake of reference, Fig 
presents the distribution of all spectra in the log T c g — log g 
plane. 



3.1. Kurucz atmospheres 

Earlier Padova isochrones were based on the Kurucz 
(1993) libraries of ATLAS9 synthetic atmospheres. As 
discussed in a series of papers by Castelli et al. (1997), 
Bessell et al. (1998), and Castelli (1999), these models 
are superseded by now. Firstly small discontinuities asso- 
ciated to the scheme of "approximate overshooting" ini- 
tially adopted by Kurucz have been corrected (cf. Bessell 
et al. 1998). Secondly no-overshooting models have been 
demonstrated to produce T e ff-colour relations in better 
agreement with empirical ones, at least for stars hotter 
than the Sun (Castelli et al. 1997). 

3.1.1. The more recent models 

In the present work, we adopt the ATLAS9 no-overshoot 
models that have been calculated by Castelli et al. 
(1997). They correspond to the "NOVER" files avail- 
able at http://cfaku5.harvard.edu/grids.html. The 
metallicities cover the values [M/H] = —2.5, —2.0, —1.5, 



4.5 



M 
O 



3.5 




□ □□□□□ 



□ □□□□□ 



□ □□□□□ 



_L 



_L 



6 4 2 

log g 

Fig. 1. Distribution of the [M/H] = spectra incorpo- 
rated in our stellar library (large symbols) in the log T c g — 
log 3 plane, compared to the position of stellar models of 
solar metallicity (small dots; these models are isochrones 
to be discussed later in Sect. 5.1). The spectra are taken 
from Castelli et al. (1997; crosses), Fluks et al. (1994; cir- 
cles), Allard et al. (2000a; squares), and pure blackbody 
(triangles). Fluks et al. spectra have been arbitrarily lo- 
cated at log g = 0. Similar distributions hold for all metal- 
licities between [M/H] = —2.5 and +0.5 (see text). 



— 1.0, —0.5, 0.0, and +0.5, with solar-scaled abundance 
ratios. A microturbulent velocity £ = 2kms _1 , and a mix- 
ing length parameter a = 1.25, are adopted. Notice that 
these models are now being extended so as to include also 
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a-enhanced chemical mixtures, which represents a poten- 
tially important improvement for our future works. 

Kurucz models cover quite well the region of the 
logics vs. log g plane actually occupied by stars, at least 
in the 3500 K < T < 50 000 K, < log 5 < 5 intervals (see 
Fig. [I]). However, it has to be extended to both lower and 
higher T c ffS, as will be detailed below. 

It is important to recall that Kurucz (ATLAS9) spec- 
tra are widely used in the field of synthetic photome- 
try, mainly because of their wide coverage of stellar pa- 
rameters and easy availability. Moreover, there are also 
good indications in the literature that these spectra do a 
good job in synthetic photometry, provided that we are 
dealing with broad-band systems. Compelling examples 
of this can be found in Bessell et al. (1998), who com- 
pares the UBVRIJHKL results obtained from the re- 
cent ATLAS9 spectra to empirical relations derived with 
the infrared flux method, lunar occultations, interferome- 
try, and eclipsing binaries. Their results indicate that the 
1998 ATLAS9 models are well suited to synthetic photom- 
etry, but for small errors, generally lower than 0.1 mag in 
colours, that we do not consider as critical. In fact, we 
are more interested in the overall dependencies of colours 
and magnitudes with stellar parameters - probably well 
represented by present synthetic spectra - than on details 
of this order of magnitude. 

Additionally, Worthey (1994) presented extensive com- 
parisons between Kurucz (1993) spectra and stars in the 
low-resolution spectral library by Gunn & Stryker (1983), 
obtaining generally a good match for wavelengths redder 
than the B pass-band. Worthey's figure 9 also presents 
a comparison between Kurucz (1993) solar spectra and 
Neckel & Labs (1984) data, with excellent results (errors 
lower than 0.1 mag) all the way from the UV up to the 
near-IR. Since the ATLAS9 1998 spectra differ just little 
from the Kurucz (1993) version (a few percent in extreme 
cases), these results are to be considered still valid. 

3.1.2. Some caveats on ATLAS9 spectra 

The previously mentioned works point to a reasonably 
good agreement between ATLAS9 spectra and those of 
real stars of near-solar metallicity, especially in the visual 
and near-infrared pass-bands. However, there are many 
known inadequacies in these spectra, which should be kept 
in mind as well. Here, we give just a brief list of the poten- 
tial problems, concentrating on those which may be more 
affecting our synthetic colours. 

ATLAS9 spectra are based on ID static and plan- 
parallel LTE model atmospheres, which use a huge 
database of atomic line data (Kurucz 1995). The line list 
is known not to be accurate: In fact, Bell et al. (1994) 
show that the solar spectra calculated using Kurucz list 
of atomic data present many unobserved lines; moreover, 
the number of lines which are too strong exceeds those 
which are too weak. The problem can be appreciated by 
looking at the high-resolution spectral plots presented by 



Bell et al. (1994), but could hardly be noticeable in low- 
resolution plots (such as in the comparisons presented in 
Worthey's 1994 figure 9, and in Castelli et al. 1997 figure 
2). 

Also, Bell et al. (2001) show that a motivated in- 
crease in the Fel bound-free opacity cause a significant 
improvement in the fitting of the solar spectrum in the 
3000 — 4000 A wavelength region, affecting the entire UV 
region as well. Such increased sources of continuous opac- 
ity are still missing in ATLAS9 atmospheres^. 

These results indicate that ATLAS9 spectra will pro- 
duce worse results when applied to (i) narrow-band pho- 
tometric systems, in which individual metallic lines can 
more significantly affect the colours, and (ii) in the UV 
region, especially shortward of 2720 A (see Bell et al. 
2001). In both cases, the errors caused by wrong atom- 
ica data are such that we can expect not only systematic 
and T e ff-dependent offsets in synthetic colours, but also a 
somewhat wrong dependence on metallicity. Clearly, these 
points are worth being properly investigated by means of 
detailed spectral comparisons. 

Regarding the present work, the above-mentioned 
problems (i) critically determine the inadequacy of 
synthetic colours computed for the Stromgren system 
(Girardi et al., in preparation), and (ii) may possibly 
cause significant errors in our synthetic HST/WFPC2 UV 
colours. 

Other potential problems worth of mention are: 

— The deficiencies in ID atmosphere models and in the 
mixing length description of convection, as compared 
to exploratory 3D model atmospheres (e.g. Asplund ct 
al. 2000). They should affect mainly the intermediate- 
strong lines with significant non-thermal broadening 
and in stars cooler than ~ 6000 K. The consequences 
of spatial and thermal inhomogeneities in the atmo- 
sphere of the Sun and Procyon have been investigated 
by Allende Prieto et al. (2001, 2002). 

— Inadequacies in the hydrogen line computations. 
Especially affected are high members of the Balmer 
series in GK stars - due to inhomogeneities in real at- 
mospheres and inadequate treatment of hydrogen line 
broadening in cool stars (see Barklem et al. 2000a, b) - 
and the core of Balmer lines and the region around the 
Paschen discontinuity and longward in A stars (Colina 
et al. 1996). In synthetic photometry, these inadequa- 
cies are expected to slightly affect the U and B mag- 
nitudes for A-F stars, but to have a greater effect in 
Stromgren indices. 

— The inadequacy of using scaled-solar metal ratios in 
giants where the abundance of CNO elements has 
been altered by dredge-up events. Similarly, the inad- 
equancy of using scaled-solar abundances for halo and 
bulge stars, instead of a-enhanced ones. 

2 In the original Kurucz (1993) spectra, a modest reduction 
of the continuous UV flux results from the adopted "approxi- 
mate overshooting" scheme. 
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Obviously, some improvement upon these points is ex- 
pected in future releases of ATLAS spectra (see Castelli 
& Kurucz 2001), and of other extended spectral grids as 
well. Fortunately, the work by Bessell et al. (1998) gives 
us some confidence that present broad-band magnitudes 
and colours (from U to K) are modelled with an accuracy 
that is already acceptable for many applications. 

Finally, we remark that some authors (Lejeune et al. 
1997, 1998) propose the application of a posteriori trans- 
formations to Kurucz (1993) spectra, as a function of 
wavelength and T g, such as to reduce the errors of the de- 
rived synthetic UBVRIJHKL photometry. In our opin- 
ion, such transformations are questionable because they 
do not correct the cause of the discrepancies - majorly 
identifiable in the imperfect modelling of absorption lines 
- and the case for applying them to stars of all surface 
metallicities and gravities is far from compelling. 

3.2. Extension to higher temperatures 

For T e g > 50 000 K, we simply assume black-body spectra. 
This is probably a good approximation for wavelengths 
A > 912 A. In fact, we find always a reasonably smooth 
transition in the computed BCs x s as we cross the T e ff = 
50 000 K temperature boundary. 

3.3. Extension to M giants 

Synthetic spectra for M giants have still many problems 
mainly in their ultraviolet-blue region - that partially 
derive from incomplete opacity lists of molecules such as 
TiO, VO and H 2 (see e.g. Plez 1999; Alvarez & Plez 
1998; Alvarez et al. 2000; and Houdashelt et al. 2000a,b 
to appreciate the state of the art in the field). 

Therefore, we prefer to use the empirical M giant spec- 
tra from Fluks et al. (1994; or "intrinsic" spectra as re- 
ferred in their paper) . They cover the wavelength interval 
from 3800 A to 9000 A. Outside this interval, the empirical 
spectra have been extendend with the "best fit" synthetic 
spectra computed by the same authors. 

However, the whole procedure reveals a problem: if we 
simply merge empirical and synthetic spectra from Fluks 
et al. (1994), the resulting synthetic B-V and U-B 
colours just badly correlate with the measured colours 
for the same stars (which were also obtained by Fluks et 
al. 1994). This problem probably derives from a bad flux 
calibration at the blue extremity of the observed spec- 
tra and/or from the imperfect match between synthetic 
and observed spectra at 3800 A. In order to circumvent 
(at least partially) the problem, we simply multiply each 
M-giant spectrum blueward of 4000 A (with a smooth 
transition in the range from 4000 A to 4800 A) by a con- 
stant, typically between 0.8 and 1.2, so that the synthetic 
colours recover the observed behaviour of the B — V vs. 
V— K data. The first two panels of Fig. || show the re- 
sults. 



Actually, Fig. g presents six different colour vs. V — K 
diagrams that are useful to understand the situation for gi- 
ants. Care has been taken in expressing data and models in 
the same photometric system, the "Bessell" UBVRIJHK 
one, that we will detail later in Sect. 4.1. For M giants, 
the empirical photometric data from Fluks et al. (1994; 
small dots) can be compared with the results of our syn- 
thetic photometry^]. Noteworthy, there is a reasonably 
good match between the synthetic and observed relations 
for most colours. This has been imposed for U — B and 
B — V, whereas is a natural result for all colours involving 
wavelengths longer than ~ 4800 A. The only clear excep- 
tion is the V — R colour, for which differences of ~ 0.4 mag 
are found for all giants of spectral type later than M4 
(V—K ib 5). The reason for this discrepancy is not clear, 
but may lie in the use of R filters with different trans- 
mission curves. Also the predictions for J — K do not fit 
well all the photometric data, somewhat failing for the 
spectral types later than M7 (V—K il 8). However, since 
these latters are quite rare, such mismatch does not pose 
a serious problem. 

For the sake of comparison, Fig. 2 also presents the 
relations obtained by means of the M-giant models from 
Houdashelt et al. (2000a), in the case of solar metallic- 
ity Together with other recent examples (e.g. Plez 1999; 
Alvarez et al. 2000), they represent state-of-the-art com- 
putations of cool oxygen- rich stellar atmospheres. As can 
be appreciated in the figure, Houdashelt et al. models re- 
produce well the empirical data as far as V—K Si 6 (spec- 
tral types earlier than M5), but start departing from these 
for cooler stars. A similar situation holds if we look at dif- 
ferent T e ff-colour relations, as can be seen in figures 13 
and 14 of Houdashelt et al. (2000a), where they compare 
their T e ff-colour relations with those obtained with Fluks 
et al. (1994) spectra and data for field giants. Also in this 
case, it seems that Fluks et al. (1994) spectra do better 
reproduce the empirical relations for the spectral types 
later than M4. 

Once we have defined the library of M-giant spectra, 
we associate effective temperatures to them by using the 
scale favoured by Fluks et al. (1994). In this scale, M giants 
cover the temperature interval from 3 850 K (MK type 
M0) to 2 500 K (MK type M10). We recall that Fluks et 
al. (1994) T e ff values are derived from a careful fitting of 
the observed spectra with synthetic model atmospheres of 
solar metallicity. Their scale is also in excellent agreement 
with the empirical one from Ridgway et al. (1980), which 
covers spectral types earlier than M6. 

After the proper T e g is attributed, each one of our 
modified spectra is completely re-scaled by a constant, so 



bol 



is 



that the total flux vs. T c g relation - 
recovered. 

Finally, we face the problem of defining the transition 
between the M-giant spectra, and the ATLAS9 ones which 



3 The JHK colours from Fluks et al. are in the ESO sys- 
tem and have been converted to the Bessell one by using the 
relations found in Bessell & Brett (1988). 
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Fig. 2. Colour vs. V — K relations for giants. The connected open circles represent the relation obtained from [M/H] = 
ATLAS9 spectra located along the T e s = 3250 + 500 logg line (the typical location for RGB stars) in the diagram 
of Fig. |]. The connected open squares correspond to the relation obtained from the M-giant spectra from Fluks et 
al. (1994), completed and modified at A < 4800 A as detailed in the text. The dashed lines represent the empirical 
relations for F0-K5 solar-metallicity giants from Alonso et al. (1999b), whereas the dotted lines are the synthetic 
relations for K0-M7 giants from Houdashelt et al. 2000a). The empirical data for M giants (small dots) are from Fluks 
et al. (1994). As far as possible, all observations have been converted to the same photometric system as used in our 
synthetic photometry (i.e. the "Bessell" UBVRIJHK system; see text). 

are available for temperatures higher than 3 500 K. To this represents fairly well the locus of low-mass red giants 
aim, it is helpful to examine Fig. 0, where we also include: in a T e g vs. \ogg diagram (see Fig. [[]). 

— the synthetic photometry for a sequence of ATLAS9 
spectra for M giants (empty circles). These spectra are 
located along the line T e g = 3250 + 500 log<?, which 
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— the mean empirical relations for F0-K5 giants of 
[Fe/H] = as derived from Alonso et al. (1999b) fitting 
formulas (dashed lines)^. 

An important fact to be noticed is that our synthetic pho- 
tometry reproduces Alonso et al. (1999b) relations in all 
colours remarkably well. 

From inspecting this and other similar plots, we can 
conclude that the mismatch between Kurucz ATLAS9 and 
Fluks et al. (1994) spectra starts at about T g = 3 850 
and increases slowly as the temperature decreases down 
to 3 500 K (i.e. from V-K~ 3.5 to V-K ~ 4.7). Hence, 
we adopt a smooth transition between these two spectral 
sources over this temperature interval. The same M giant 
spectra are assumed for all metallicities. 

The complete procedure ensures reasonable colour vs. 
V—K relations for all giants of near-solar metallicity 
(Fig. ||). Nevertheless, this kind of approach cannot be 
completely satisfactory, first because the original Fluks et 
al. (1994) spectra have been artificially corrected at wave- 
lengths shorter than 4800 A in order to produce reasonable 
B — V and U — B, and second because we do not dispose 
of similar M-giant spectra for metallicities very different 
from solar. Better empirical and theoretical spectra for M 
giants seem to be urgently needed. Anyway, in the context 
of the present work the problem is not dramatic because 
M giants cooler than T e g ~ 3 500 K are only found in the 
RGB-tip and TP-AGB phases of high metallicity stellar 
populations, and constitute just a tiny fraction of the num- 
ber of red giants. The problem could be critical, instead, 
when we consider integrated properties of stellar popula- 
tions, because M giants, despite their small numbers, have 
high luminosities and contribute a sizeable fraction of the 
integrated light. 



better opacity lists and dust formation. Dust can signifi- 
cantly affect the coolest atmospheres, corresponding to 
dwarfs of spectral types L and T. 

The selected spectra cover the T e g intervals: 

- from 4 000 K to 2 800 K ("AMES" models) for metal- 
licities [M/H] = 0.0, -0.5, -1.0, -1.5, and both 
log g = 5.0 and 5.5; 

- from 2 800 K down to 500 K ( ' AMES-dusty" models) 
only for [M/H] = 0.0, and log g values between 3.5 and 
6.0. 

These spectra are presented with a extremely high reso- 
lution, that by far exceeds the one necessary in our work. 
Thus, we have convolved the flux per unit frequency F v 
with a Gaussian filter of er„ = 2.4 x 10~ 18 Hz, that corre- 
sponds to a FWHM of 20 A at A = 5550 A. The resulting 
spectra were then reported to the same grid of wavelengths 
of Kurucz' spectra. 

We find that there is a good agreement between 
ATLAS9 and BDdustyl999 spectra in the T off range be- 
tween ~ 3 800 K and 4 000 K. Then, we set the transition 
between ATLAS9 and BDdustyl999 spectra at ~ 3 900 K. 
This choice guarantees smooth T c g vs. colour relations for 
dwarfs. 

4. Photometric systems available 

In this section we present the basic information regarding 
the filter transmission curves and zero-points for each pho- 
tometric system. As a reference to the discussion, Fig. || 
presents the filter sets under consideration, as compared 
to the spectra of a hot (Vega), an intermediate (the Sun), 
and a cool star (an M5 giant). 



3.4. Extension to M+L+T dwarfs 



4.1. Johnson-Cousins-Glass system 



Although the modelling of cool dwarfs atmospheres 
presents challenges comparable to those found in late-M 
giants (e.g. the inadequacy of TiO and H 2 line lists, and 
dust formation; see Tsuji et al. 1996, 1999; Leggett et al. 
2000), present results compare reasonably well with obser- 
vational spectral data (see e.g. figure 9 in both Leggett et 
al. 2000 and 2001). A review on the subject can be found 
in Allard et al. (1997). 

An extended library of synthetic spec- 
tra for cool dwarfs (of types M and later) 
is provided by Allard et al. (2000a; see 
f tp : //f tp . ens-Lyon . f r/ pub/users/CRAL/f allard) . 
We use their set of "BDdustyl999" atmospheres (see also 
Chabrier et al. 2000; Allard et al. 2000b, 2001), that 
should supersede the "NextGen" models from the same 
group (Hauschildt et al. 1999) due to the consideration of 



4 Alonso et al. (1999b) V—I and V — R colours have been 
transformed from Johnson to Johnson- Cousins systems using 
relations from Bessell (1979), whereas infrared colours were 
converted from TCS to "Bessell" systems using the relations 
from Alonso et al. (1998). 



Aiming to reproduce the Johnson-Cousins-Glass system, 
we adopt the filter pass-bands indicated by Bessell (1990; 
for Johnson-Cousins UBVRI) and Bessell & Brett (1988; 
for JHK) . We also apply their prescription for computing 
the U—B colour by means of a slightly modified pass-band 
BXqo, instead of the normal B one. Moreover, in order 
to better recover the original system, we adopt energy 
instead of photon count integrations (see Sect. 2.1). 

It is worth recalling that these pass-bands repre- 
sent just one specific version of the "standard" Johnson- 
Cousins-Glass system, that may differ from filter systems 
in usage at several observatories. The Bessell & Brett 
(1988) JHKLM pass-bands, for instance, represent an 
effort in the direction of homogenizing several different 
near-infrared systems (SAAO, ESO, CIT/CTIO, MSO, 
AAO, and Arizona). In Bessell (1990) and Bessell & Brett 
(1988), the reader can find a set of useful fitting relations 
between colours in the several original systems. 

As previously mentioned, Johnson-Cousins-Glass is es- 
sentially a VEGAmag system. We fix the zero-points by 
assuming that Vega has apparent magnitudes equal to 0.03 
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Fig. 3. The filter sets used in the present work. From top to bottom, we show the filter+detector transmission curves 
S x for the systems: (1) HST/NICMOS, (2) HST/WFPC2, (3) Washington, (4) ESO/EMMI, (5) ESO/WFI UBVRIZ 
+ ESO/SOFI JHK, and (6) Johnson-Cousins-Glass. All references are given in Sect. ||. To allow a good visualisation 
of the filter curves, they have been re-normalized to their maximum value of S\. For the sake of comparison, the 
bottom panel presents the spectra of Vega (AOV), the Sun (G2V), and a M5 giant, in arbitrary scales of F\. 



in all UBVRIJHK bands, i.e. we impose all colours to 
be null. Notice that our definition is very similar to the 
Bessell et al. (1998) one, who adopt Vega colours differing 
from zero by just some thousandths of a magnitude (see 
their table Al). 



4.2. Instruments on board of HST 

In the context of the present work, the distinctive feature 
of HST photometry is the use of ABmag and STmag sys- 
tems, which greatly simplifies the definition of zero-points. 
WFPC2 and NICMOS observations can also be expressed 
in VEGAmag magnitudes. 
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4.2.1. WFPC2 

As for the WFPC2, we produce bolometric corrections 
and magnitudes in the F170W, F218W, F255W, F300W, 
F336W, F439W, F450W, F555W, F606W, F702W, 
F814W, and F850LP filters. Similar tables have been pro- 
duced in previous papers (Chiosi et al. 1997; Salasnich et 
al. 2000); the present ones differ just in minor details. 

The transformations have been computed in STmag, 
ABmag, and VEGAmag systems. Whereas the STmag 
and ABmag systems can be straightforwardly simulated, 
the VEGAmag one deserves some comments. 

According to the SYNPHOT package distributed with 
the STSDAS software, in a VEGAmag system Vega should 
have apparent magnitudes equal to zero in all pass-bands. 
However, the most widely used calibration of WFPC2 pho- 
tometry comes from Holtzman et al. (1995), who adopt a 
set of zero-points in which Vega has apparent magnitudes 
U = 0.02, B = 0.02, V = 0.03, R = 0.039, I = 0.035. 
Accordingly, we choose this latter definition, and impose 
that for the WFPC2 filters that correspond to UBVRI in 
wavelength, the synthetic magnitudes of Vega have these 
same values. For filters of intermediate wavelength, we 
adopt a linear interpolation between these values, whereas 
for bluer and redder filters a Vega magnitude equal to zero 
is assumed. 

As for the S\ functions, we use the pre-launch pass- 
bands kindly provided by Jon Holtzman (sec Holtzman et 
al. 1995). We recall that, owing to the presence of contam- 
inants inside the WFPC2 (see Baggett & Gonzaga 1998; 
Holtzman et al. 1995), S\ changes slowly with time, espe- 
cially for UV filters. To cope with this, observers usually 
apply small corrections to the definition of the instrumen- 
tal magnitudes, in order to bring the magnitudes back to 
the original conditions (see Holtzman et al. 1995). This 
justifies our use of pre-launch pass-bands instead of the 
present-day ones provided by SYNPHOT. 

4.2.2. NICMOS 

For NICMOS filters, we compute bolometric corrections 
and absolute magnitudes in the ABmag, STmag, and 
VEGAmag systems. For the moment, calculations are 
limited to the three most frequently used filters, i.e. 
F110W, F160W, and F205W. The pass-bands come from 
SYNPHOT, and have been kindly provided by Don Figer. 

NICMOS observations are also frequently expressed in 
units of milli-Jansky (mJy). The conversion between AB 
magnitudes and mJy is straightforward: 

mAB,^(mag) = -2.51og/ l/ (mJy) - 16.4 . (16) 

4.3. ESO Imaging Survey (EIS) filter sets 

The EIS survey (Renzini & da Costa 1997; da Costa 2000) 
aims at providing a large database of deep photometric 
data, among which the astronomical community could se- 
lect interesting targets for VLT spectroscopy. The survey 
is conducted at several different instruments at ESO: 



4.3.1. WFI 

The Wide Field Imager (WFI) at the MPG/ESO 2.2m 
La Silla telescope provides imaging of excellent quality 
over a 34' x 33' field of view. It contains a peculiar set 
of broad-band filters, very different from the "standard" 
Johnson-Cousins ones. This can be appreciated in Fig. || 
notice in particular the particular shapes of the WFI B 
and / filters. Moreover, EIS makes use of the WFI Z filter 
which does not have a correspondency in the Johnson- 
Cousins system. 

Given the very unusual set of filters, the importance of 
computing isochrones specific for WFI is evident. This has 
benn done so for the broad WFI filters U (ESO#841), B 
(ESO#842), V (ESO#843), R (ESO#844), I (ESO#845), 
and Z (ESO#846), that - here and in Fig. || - are referred 
to as UBVRIZ for short. 

Bolometric corrections have been computed in the 
VEGAmag system assuming all Vega apparent magni- 
tudes to be 0.03, and in the ABmag system, which is 
adopted by the EIS group. The photometric calibration 
of EIS data is discussed in Arnouts et al. (2001). 

It is very important to notice that any photometric 
observation performed with WFI that makes use of stan- 
dard stars (e.g. Landolt 1992) to convert WFI instrumen- 
tal magnitudes to the standard Johnson-Cousins UBVRI 
system, will not be in the WFI VEGAmag system we are 
dealing with here. Instead, in that case, we should bet- 
ter compare the observations with the isochrones in the 
standard Johnson-Cousins system, or, alternatively, ap- 
ply colour transformations to our WFI isochrones so that 
they reproduce the T e ff-colour sequences of dwarfs and gi- 
ants in the Johnson-Cousins system. This aspect will be 
better illustrated Sect. |[ 

4.3.2. SOFI 

SOFI is a near-infrared imager and grism spectrograph 
at the ESO NTT 3.6m telescope. It is used to comple- 
ment WFI observations in near-infrared pass-bands. More 
specifically, it disposes of J, H and K s filters that have 
some similarity to Glass JHK (see Fig. ||). 

VEGAmag magnitudes are computed assuming the 
same as before, i.e. that Vega apparent magnitudes are 
0.03. 

4.3.3. EMMI 

EMMI is a visual imager and grism spectrograph also at 
the NTT. The EIS project uses a set of its wide filters, i.e. 
B w , I w , R w , and V w , together with the R filter. The pass- 
bands have been kindly provided by S. Arnouts. As before, 
VEGAmag magnitudes are computed assuming that Vega 
has 0.03 mag in all filters. 
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4.4. Washington system 



5.1. An extended set of isochrones 



The Washington system, originally defined by G. 
Wallerstein and developed by Canterna (1976), has been 
more and more used since Geisler (1996) defined a set of 
CCD standard fields. In the present paper, we reproduce 
the Washington CCD system (filters C, M, T x and T 2 ) 
by adopting the transmission curves as revised by Bessell 
(2001). 

Following Geisler (1996), the Washington T\ filter has 
been frequently replaced by the Kron-Cousins filter i?, 
which presents a similar transmission curve and is avail- 
able in almost all observatories. Additionally, Holtzman 
et al. (in preparation) suggest the use of BVI colours 
together with Washington ones for breaking the age- 
metallicity degeneracy of stellar populations in colour- 
colour diagrams. For these reasons, in all tables that deal 
with the Washington CMT1T2 system, we also insert the 
information for the Johnson-Cousins BVRI filters. 

Finally, following Geisler (private communication), the 
zero-points should be well represented by a VEGAmag 
system. We assume Vega has magnitude 0.03 in all filters. 

4.5. Other planned systems 

Forthcoming papers of this series will be dedicated to 

— other traditional systems like Stromgren, Thuan- 
Gunn, Vilnius, etc.; 

— systems corresponding to some successful and re- 
cent/ongoing observational campaigns, like Hipparcos, 
MACHO, DENIS, 2MASS, and SDSS; 

— newly-proposed systems, like the ones for the future 
astrometric mission GAIA. 

Since the procedure for computing BCs x tables is rela- 
tively easy, provided the necessary information - i.e. filter 
transmission curves and a zero-point definiti on t hat corre- 
sponds to one of the cases discussed in Sect. 2J3 - is given, 
we can provide isochrone tables for any photometric sys- 
tems upon request. 

5. Available stellar tracks and isochrones 

The tables of bolometric corrections here described have 
been primarily constructed to be applied to the Padova 
database of stellar evolutionary tracks and isochrones^. 
These latter have been described in several previous pa- 
pers, and a complete description of them is beyond the 
scope of this work. 

In the following, our intention is just to briefly men- 
tion the sets of isochrones which are the most useful for 
comparisons with observed photometric data, and to first 
mention some data that has not been published in prece- 
dence. A summary table of the available material is pre- 
sented in Tabic |[ 

5 Of course, they can also be applied to any other set of 
stellar data in the literature. 



Girardi et al. (2000) computed a set of low- and 
intermediate-mass stellar tracks which supersedes those 
previously used in Bertelli et al. (1994) isochrones. 
Additional models (unpublished) have been recently com- 
puted for initial chemical composition [Z = 0.0001, Y = 
0.23]. Thus, we have created a new set of isochrones com- 
bining the latest low- and intermediate mass tracks (in the 
range 0.15 — 7.0 Mq) with the formerly available massive 
ones. The full references are given in Table ^. 

In all cases in which "2000" and "1994" tracks have 
been combined, we find excellent agreement between the 
relevant quantities (lifetimes, tracks in the HR diagram) at 
the transition mass of 7 — 8 Mq . This is explained consid- 
ering that the intermediate- and high-mass models share 
the same prescription for convection, and have interior 
opacities dominated by electron scattering (which has not 
been changed in the meanwhile). For the solar metallic- 
ity, we have combined tracks presenting slightly different 
initial metallicities - [Z = 0.019, Y = 0.273] in Girardi et 
al. (2000), and [Z = 0.020, Y = 0.280] in Bressan et al. 
(1993) - without finding any significant discontinuity. 

Finally, to this extended set we have included the 
Marigo et al. (2001) isochrones for zcro-metallicity stars. 

5.2. Overshooting vs. classic models 

For solar metallicity and in the mass range 0.15 — 
7.0 Mq, Girardi et al. (2000) presents an additional 
set of tracks and isochrones computed with the classical 
semi-convective prescription for convective borders. Then, 
the two [Z = 0.019, Y = 0.273] sets are useful if one 
chooses to compare classical with overshooting models. 
The same kind of work is being extended to metallici- 
ties [Z = 0.004, Y = 0.240] and [Z = 0.008, Y = 0.250] 
(Barmina et al. 2002), and will be included in the database 
when completed. 

5.3. Solar-scaled vs. a-enhanced models 

Salasnich et al. (2000) presents new models for 4 dif- 
ferent metallicities ([V = 0.250, Z = 0.008], [Y = 
0.273, Z = 0.019], [Y = 0.320, Z = 0.040] and [Y = 
0.390, Z = 0.070]), computed both with scaled-solar and 
alpha-enhanced distributions of metals. The tracks cover 
the mass range from 0.15 to 20 Mq. They represe nt a 
valid alternative to the isochrones referred to in Sect. |5.l| , 
but do not cover the low-metallicity interval. Anyway, it 
has been demonstrated by Salaris et al. (1993; see also 
Salaris & Weiss 1998; and VandenBerg 2000), that for low 
metallicities one may safely use scaled-solar stellar models 
instead of a-enhanced ones of same Z. 



5.4. Simple vs. more detailed TP-AGB models 

All the tracks and isochrone sets above mentioned, include 
the complete TP-AGB phase as computed with a simple 
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Table 1. Stellar tracks used in the different set of isochrones. 
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1 S indicates a solar-scaled distribution of metals, A indicates an a-enhanced one. The adopted metal abundance ratios 
are specified in Salasnich et al. (2000). 

2 References for tracks and isochrones: Be94 = Bertelli et al. (1994, A&AS 106, 275); Br93 = Bressan et al. (1993, A&AS 100, 
647); Fa94a = Fagotto et al. (1994a, A&AS 104, 365); Fa94b = Fagotto et al. (1994b, A&AS 105, 29); Gi96 = Girardi et al. 
(1996, A&AS 117, 113); GiOO = Girardi et al. (2000, A&AS 141, 371); GiOl = Girardi (2001, unpublished); MaOl = Marigo et 
al. (2001, A&A 371, 152); MG01 = Marigo & Girardi (2001, A&A 377, 132); SaOO = Salasnich et al. (2000, A&A 361, 1023). 

3 G means a simple synthetic evolution as in Girardi & Bertelli (1998), whereas M stands for more detailed calculations as in 
Marigo (2001; and references therein). 

4 O means a model with overshooting (see GiOO for all references), whereas C corresponds to classical semi-convective models. 



synthetic algorithm (Girardi & Bertelli 1998). In paral- 
lel, Marigo (2001, and references therein) has developed a 
much more sophisticated code for synthetic TP-AGB evo- 
lution, which includes crucial processes such as the third 
dredge-up and hot-bottom burning. Sets of complete TP- 
AGB tracks have been so far presented for metallicities 
[Y = 0.240, Z = 0.004], [Y = 0.250, Z = 0.008], and 
[Y = 0.273, Z = 0.019] (Marigo 2001). A set of isochrones 
has been generated by combining these detailed TP-AGB 
tracks with the previous evolution from Girardi et al. 
(2000); they are presented in the appendix of Marigo & 
Girardi (2001). 

These isochrones represent a u seful alternative to the 
isochrones referred to in Sect. Bd, any time the TP-AGB 
population is under scrutiny - for instance, when we have 
near-infrared photometry of objects above the RGB-tip. 
Marigo TP-AGB models are being extended to cover a 
larger metallicity interval. 



6. Discussion and concluding remarks 



6.1. Summary 

This work is dedicated to the presentation of theoretical 
isochrones in several photometric systems. It represents 
the continuation of a wide project of the Padova group, 
started with Bertelli et al. (1994) and then carried on by 
Chiosi et al. (1997), and more recently by Salasnich et al. 
(2000). It starts describing the formalism for converting 
synthetic stellar spectra into tables of bolometric correc- 
tions (Sect. |2|), in such a way that it can be easily applied 
to different photometric systems, and with the possibility 
of including extinction in a self-consistent way. 

Then, we describe the assemblage of an updated li- 
brary of stellar spectra (Sect. ||). The library is quite 
extended in T c ff and log g, and includes the crucial de- 
pendence of spectral features on metallicity. Of course, 
it suffers from some limitations: very hot stars and M- 
giants are not included among synthetic spectra, a situ- 
ation which can be remediated by using blackbody and 
empirical spectra, respectively. The strinkingly different 
spectra of carbon stars will also have to be considered in 
the future (Marigo et al., in preparation). These are prob- 
ably the points where the models can be most improved. 
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Fig. 4. Comparison of Girardi et al. (2000) isochrones transformed to Johnson-Cousins-Glass magnitudes and colours 
by using either present relations (continuous lines) or Bertelli et al. (1994) ones (dashed lines), in several CMDs. The 
isochrones have solar metallicity (Z — 0.019) and ages 10 8 , 10 9 , and 10 10 yr (from top to bottom). 



Moreover, several problems may be affecting the ATLAS 9 
synthetic spectra we are using t o sim ulate the broad-band 
colours of most stars (see Sect. Although such spec- 



tra have been demonstrated to be suitable for synthetic 
photometry (mainly for the Johnson-Cousins-Glass sys- 
tem; e.g. Bessell et al. 1998), their accuracy has still to 
be sistematically evaluated for stars of all metallicities, 
temperatures and gravities. Anyway, we assume that they 
are good enough for simulating broad-band photometric 
systems in the visual-infrared wavelength region, whereas 



expect that the results in narrow-band systems, and in the 
ultraviolet pass-bands, will be affected by more significant 
errors. Another important aspect of synthetic spectra is 
that they are usually computed for scaled-solar chemical 
compositions, whereas the extension to peculiar and a- 
enhanced mixtures would be of high interest. This latter 
problem will probably be alleviated in a near future, with 
the release of extensions to ATLAS9 spectra. 

From the spectral library, we derive bolometric correc- 
tions for each pass-band mentioned in Sect. |J, and apply 
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Fig. 5. Comparison of Girardi et al. (2000) isochrones as transformed to the WFPC2 VEGAmag system using either 
present relations (continuous lines) or Salasnich et al. (2000) ones (dashed lines). The isochrone ages and metallicities 
are the same as in Fig. f|. 



them to the Padova isochrones. In practice, only in few 
cases we present newly-constructed isochrones: the bulk 
of isochrone data is already described in previous papers 
by our group, and it is only the transformation from the- 
oretical quantities to absolute magnitudes that changes 
compared to past releases. Sect. summarizes the basic 
characteristics of the different sets of isochrones, indicat- 
ing the full references. 



6.2. Comparison with previous works 

It is important to illustrate the differences between 
the present and previous transformations. The previous 
ones for UBVRIJHK are fully described in Bertelli et 
al. (1994), and were adopted by Girardi et al. (2000), 
Salasnich et al. (2000) and Marigo et al. (2001); for 
HST/WFPC2 photometry, they are the ones described in 
Salasnich et al. (2000). 

The situation for Johnson-Cousins-Glass is tentatively 
illustrated in Fig [|, which compares a set of Girardi et al. 
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(2000) isochrones, transformed according to both present 
(continuous lines) and Bertelli et al. (1994; dashed lines) 
transformations. We point out that: 

1. For stars hotter than T e g ~ 4000 K, present trans- 
formations are very similar to those of Bertelli et al. 
(1994). The differences amount to just a few hun- 
dredths of a magnitude over most regions of the 
CMD, including the entire main sequence and sub- 
giant branch, and most of the RGB. They can be en- 
tirely attributed to the slightly different pass-bands 
and zero-points, and to the use of more recent ATLAS9 
"NOVER" atmospheres instead of Kurucz (1993) ones. 

2. A somewhat similar situation holds for dwarfs cooler 
than T ff ~ 4 000 K (see bottom end of the main se- 
quence in all panels, for My is 7, and Mk <; 5). The 
differences are generally small and can be attributed to 
the change from Kurucz (1993) to Allard et al. (2000a) 
spectra. The exceptions are U — B and J—H colours, 
for which the differences between the two versions be- 
come sizeable. 

3. For giants cooler than 4 000 K, present transforma- 
tions become very different. This can be noticed in 
the upper-right corner of all diagrams. M-giants corre- 
sponding to the RGB-tip and TP-AGB, are now seen 
to fade by some magnitudes in V, due to a sort of rapid 
increase of visual BCs at T e g 5s 3 500 K. This effect is 
caused by the use of Fluks et al. (1994) spectra and 
their T c g vs. V — K scale. Such a bending of the RGB 
is indeed observed in CMDs of old metal-rich clusters 
(see e.g. Ortolani et al. 1990; and Rich et al. 1998), 
and seems to be better reproduced now than with pre- 
vious transformations. Other differences appear in all 
colours: the most remarkable are the excursion of M- 
giants of latest type towards much bluer U — B, and 
the much smoother behaviour now obtained for J — H 
and H-K. 

A quite similar situation holds for HST/WFPC2 pho- 
tometry, as illustrated in Fig. ||. This time, we compare 
the same isochrones as transformed with present (contin- 
uous lines) and Salasnich et al. (2000; dashed lines) trans- 
formations. It is evident that the present transformations 
ensure a more continuous behaviour of the colours for all 
low-temperature stars (both dwarfs and giants). 

From the plots at the top row of Fig. |^, one can also 
appreciate the unusual appearance of isochrones in CMDs 
that involve ultraviolet WFPC2 pass-bands: Notice for 
instance that in F170W, F218W, F255W and F330W 
magnitudes, giants may be fainter than turn-off stars. 
Isochrones in the F218W vs. F170W-F218W and F255W 
vs. F218W— F255W diagrams are even "twisted", because 
the T e ff vs. colour relations are not monotonic for these 
filters. These effects are related to the presence of a red 
leak in the ultraviolet HST filters (for both the present 
WFPC2 and the former FOC camera), and are extensively 
discussed by Yi et al. (1995) and Chiosi et al. (1997). 

As a consequence of the great similarity between 
present and previous UBVRIJHK and HST/WFPC2 



transformations, for most colours and over a large portion 
of the HR diagram, most results derived from previous 
Padova isochrones are not expected to change. Exceptions 
may show up for works that are concerned with the 
photometry of the reddest giants, with T c g Ss 3 500 K 
(B — V <; 1.5), or that deal with low-mass main-sequence 
stars in the U — B and ultraviolet colours. 

6.3. New results 

The greatest improvement of the present database is in the 
presentation of Padova isochrones in several photometric 
systems for which they were not available so far - including 
the case of brand-new systems. Three examples of theis 
kind are given in Figs. |, @ and | 

First, Fig. ^| presents the isochrones in NICMOS 
ABmag system. In a VEGAmag system, NICMOS 
isochrones would look similar to their equivalent Johnson- 
Cousins-Glass ones, shown in Fig. || In the ABmag sys- 
tem, however, they appear shifted to quite different colour 
and magnitude intervals. 

Figure. |t| illustrates how Padova isochrones look like in 
the Ti vs. C - Tl CMD of Washington photometry, both 
for varying age at constant metallicity (left panel), and 
for varying metallicity at constant age (right one). The 
striking feature in these plots is the excellent separation 
in metallicity offered by the C — T\ colour, from the main 
sequence up to red giant phases. This feature, combined to 
the excellent throughput in the C filter (Fig. ||), is among 
the advantages that make the Washington system a very 
competitive one if compared to Johnson-Cousins (see also 
Paltoglou & Bell 1994, and Geisler & Sarajedini 1999). 

Preliminary comparisons point to a good agreement 
between our Washington isochrones and real data for LMC 
fields from Bica et al. (1998). Just to mention an example, 
we notice that C — T\ for giants "saturates" at ~ 3.4, both 
in the models and in the LMC data. 

An example of "new" photometric system is provided 
by the WFI, which has broad-band filters very differ- 
ent from Johnson-Cousins ones. To illustrate the effect in 
colours, Fig. |8| shows exactly the same isochrones as seen 
in BVI CMDs using either Johnson-Cousins or WFI fil- 
ters, and applying in both cases the VEGAmag definition 
of zero-points. The differences are striking. In particular, 
since the BV WFI filters represent a wavelength baseline 
shorter than the Johnson ones, they provide a more mod- 
est separation of stars in B — V colour. It is evident from 
this plot that the normal Johnson-Cousins isochrones can- 
not be used to interpret WFI data that has been converted 
to VEGAmag or ABmag systems, as for most of EIS data 
(e.g. Arnouts et al. 2001; Groenewegen et al. 2002)0. 



Notice, however, that part of the data released from EIS - 
namely the pre-FLAMES survey - is in fact converted into a 
standard Johnson- Cousins system (see Momany et al. 2001). 
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Fig. 6. Isochrones in the CMDs 
of NICMOS ABmag photometry. 
Ages and metallicities are the 
same as in Fig. [I| 
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Fig. 7. Isochrones in the 7\ vs. 
C — T\ plane of Washington pho- 
tometry. Left panel: from top to 
bottom, a sequence of Z = 0.008 
isochrones with ages 10 7 , 10 8 , 
10 9 , and 10 10 yr. Right panel: 
from left to right, a sequence of 
14 Gyr old isochrones with metal- 
licities Z = 0.0001, 0.0004, 0.001, 
0.004, 0.008, 0.019, and 0.030. 



6.4. Retrieval of electronic tables 

All the data here mentioned are available at the WWW 
site http://pleiadi.pd.astro.it. The database al- 
ready includes a very large number of files, and is expected 
to increase further as we publish data for other photomet- 
ric systems. Thus, it is hard to describe here both the 
structure of the database, and the content of each file. 
Moreover, this kind of information is probably useful just 
to whom actually accesses the database. Thus, we opt to 



provide all the relevant information in readme . txt files 
inserted in the database. 

To the general reader, suffice it to briefly mention the 
kind of data which is available: 

— Tables of bolometric corrections for each metallicity: 
they contain the quantities BCs x for each filter, and 
as a function of stellar T c s and log g. Metallicities avail- 
able are [Fe/H] = -2.0, -1.5, -1.0, -0.5, 0, +0.5. The 
values of T e ff and log g are not exactly the same for all 
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Fig. 8. Comparison between the same set of Z = 0.019 isochrones as seen in the UBVI CMDs using either Johnson- 
Cousins (dashed lines) or WFI filters (continuous lines). VEGAmag systems are used in both cases. Ages are the same 
as in Fig. ||. 



metallicities, but correspond quite well to the regions 
indicated in Fig. [|. 

— Tables of isochrones: they include all metallicities 
and cases indicated in Sect. |^. The various photo- 
metric systems are separated in different directories. 
For each isochrone file, the information and structure 
are the same as already presented in Girardi et al. 

(2000) , Salasnich et al. (2000), and Marigo & Girardi 

(2001) , with the obvious difference that instead of 
UBVRIJHK absolute magnitudes, in each file we 
tabulate the absolute magnitudes for the photometric 
system under consideration. 

— Tables of integrated colours of single-burst stellar 
populations: a table of this kind is present for each 
isochrone table. They provide integrated magnitudes 
in each pass-band, as a function of age. 
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